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Abstract 

We present analytic expressions for the single particle excitation energies of the 8 
quasi-particles in the lattice Eg Ising model and demonstrate that all excitations have an 
extended Brillouin zone which, depending on the excitation, ranges from < P < 4n to 
< P < 12n. These are compared with exact diagonalizations for systems through size 10 
and with the Eg fermionic representations of the characters of the critical system in order 
to study the counting statistics. 



1. Introduction 

Several years ago Bazhanov, Nienhuis and Warnaar [|J demonstrated that the the 
solution of the dilute A 3 model of Warnaar, Nienhuis and Seaton @-0 may be reduced 
to the Bethe's Ansatz equation based on the group Eg given in ||-[[|. This thus provides 
a lattice realization of the Eg continuum field theory found by Zamolodchikov J7|-[f§ in 
1989 to be in the same universality class as the Ising model in a magnetic field at T = T c . 

The reduction in [|I| to the Eg equations is based on a conjecture for the allowed string 
type solutions of the Bethe's Ansatz equations for the dilute A3 model and this conjecture 
was motivated by a numerical study of the system at criticality in the sector / = 1 which 
contains the ground state. This study gives a set of counting rules for the spectrum which 
exactly reproduces the Eg fermionic representation of the character x^i^iQ) °f the Ising 
model which was conjectured in H-[]10|] and proven in its polynomial generalization in 



However, there is as yet no satisfactory conjecture for the string content at criticality 
of all states in the sectors / = 2 and / = 3 (which correspond to the characters X^i^il) 
and Xi 3 ^ 4 ^)) and a recent study |12[-[T^] of numerical solutions to the non-critical Bethe's 



equations at zero momentum indicates that the string identification may at times change 
as the system moves away from criticality. We make a two fold study of these questions 
here by numerically computing the single particle excitations for finite chains of sizes up 
through 10 and by comparing these results with the single particle dispersion relations 
which are obtained by extending the computations of [Q. We find that all eight quasi- 
particles exhibit the phenomenon of an extended Brillouin zone such as was first explicitly 



seen in the anti-ferromagnetic three state Potts spin chain [[E_4| and that there are no 
genuine one particle states in the two sectors / = 2 and 3. 
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We present the dispersion relations in Sec. 2, the numerical study in Sec. 3 and discuss 
our findings in Sec. 4. We follow the notation of ref. and refer the reader to that paper 
for the explicit Boltzmann weights which define the model. We will here be concerned only 
with the Hamiltonian version of the model. 



2. Single Particle Dispersion Relations and Characters 



In [0 and Q the single particle excitations were studied in the limit L — > oo and in 
particular the Fourier transform of the single particle excitations in sector / = 1 is com- 
puted. We have inverted the transform to find the results given below. The computations 
involve some tedious algebra and moreover for the non-critical case the required integrals 
do not seem to be in the literature. The details will be given elsewhere. 

For the critical case we have as a function of rapidity u 



ej(u) 



n + 2 arctan 
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A dP [ — = A cosh u 
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(2.1) 



du 



-~* cosh u — cos 2 (a7r/30) 



where A = 32/15 , — oo < u < oo, the number thirty has the significance of being the 
Coxeter number of Eg and a takes on the following values as a function of j (we use the 
labeling of the Cartan matrix E$ shown in Fig. 1 and indicate the identification of particles 
of [|]J in parentheses) 
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The numbers a in this table agree with the corresponding numbers for the scattering of 
a particle of type 1 with a particle of one of the 8 types of the E 8 field theory as given 
in |y| and have group theoretical and geometrical interpretations |16|-|T^]. They also 



follow from (4.19) of 0. The numbers AP are 2ir times the elements in the first column 
of the inverse Cartan matrix. 

For the non-critical case in which the nome qs of the Boltzmann weights of U (which 
is essentially the magnetic field) is not zero we have the following generalization which 
reduces to ( |2.1|) when q — > (and the modulus k — > 0) 



v-^ [ / % sn iu \ 
Pjiu, q) = > 7r + 2 arctan — ; ; — - 



(2.3) 

e, («, q) = A(q) = A(q) dmu cn iu ^ cn2z ^ a f/ ( 1 a ^ /15) 

where q = q^^ 15 , A(q) = 2AK(q)/ir, K(q) is the complete elliptic integral of the first kind 
and — K' < u < K' . This expression is in fact a universal form for all the Bethe's Ansatz 
models based on a simply laced Lie algebra, the only difference being that the values of a, 
here given by (|2.2|) , vary from model to model. 

The minimum in the energy occurs u = ±K' where P = (or the value AP given in 
(p|)). Thus we find 

ej(P = 0,q) = kA(q)^2sn(aK/15) (2.4) 

a 

and in particular 

lime i (P = 0,q)/kA(q) = ^ sin(a7r/30). (2.5) 

a 

It may be verified that these values coincide with the components of the Perron- Frobenius 



eigenvector of the Eg Cartan matrix (as given, for example, in [|19 
When q — > 1 (fc — >• 1) we find that 



(u,q)^A(q)J2l (2-6) 



e. 



where the normalizing constant is diverging. Thus using (|2.2|) we see that the single 
particle states for particles 1 and 2 become degenerate as q — > 1 (and similarly for particles 
3 and 4). For the remaining particles degeneracies with multi-particle states also occur. 

In principle u can be eliminated between the two expressions in (|2.3|) to produce a 
polynomial relation between e j and Pj . However here we have done the elimination numer- 
ically and present the results in Fig. 2 for qs = and in Fig. 3 for qs = 0.2. In Fig. 3 we 
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see that the labeling we have used corresponds to the ordering of eigenvalues at P = 
which is the same for all q > 0. It is to be explicitly remarked that the dispersion relations 
do not have the restriction < P < 2tt but in all cases have the larger momentum range 
< P < AP where AP is given in ( |2.2| ) . This is the phenomenon of the extended Brillouin 
zone scheme. This is a very general phenomenon in integrable models first explicitly seen 
14] for the three state anti-ferromagnetic Potts spin chain. 



m 



We also wish to make contact with the characters of conformal field theory and finite 
size computations. The fermionic representation of characters of the Ising model in the Eg 
basis is given 



10 1 in terms of the fermionic form 
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where (q) n = f7" =1 (l — qi), Ce 8 is the Cartan matrix of the Lie algebra E 8 given by the 
incidence matrix of Fig. 1 (where we use the labeling of [0) and we trust that q will not 
be confused with the nome of the elliptic functions. Explicitly 
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proven in [II] that 
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(2.9) 



and in [|10| it was conjectured on the basis of computer studies that 
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(2.10) 
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We also note that the character ( |2.9|) is derived from the lattice quasiparticle spectrum 



(2.11) 



3 = 1 i=l 



8 Hj 
3 = 1 i=l 



(2.12) 



where ej(P) are single particle dispersion relations, we have the fermionic restriction 



P] ^ for i^k and all j 



(2.13) 



and the momenta are chosen from the set 
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An 



P] E \ Pp n (n), P™ n (n) + — , Pf in (n) + — , P™ x (n) 
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with 



and 



(nci), + i(i-%: 



C X (n) = -Pf in (n) + 2 7 r(C, 1 ) lj , 



(2.15) 



(2.16) 



3. Finite Size Study 

We have numerically determined the Hamiltonian eigenvalues for chains up through 
size 10. We have determined the sector I from the property that at q = eigenvectors in 
the sector I = 2(1 = 1,3) are antisymmetric (symmetric) under the interchange 1 <-> 3 of 
the states in the basis of 0, and by a similar symmetry property under the interchange 
2 <-> (1 + 3)/v2 that distinguishes / = 1 from / = 3. We present these data in Figs. 4-7 
where we also compare with the L — > oo formula ( |2.3j ). 

In Fig. 4 we plot at q = the eigenvalues which behave as 2A(q) as q — > 1. These 
include all the single particle states for particles 1 and 2 in the sector / = 1 which are 
allowed by (|2.14| ). These states are indicated with triangles. The remaining states are 
to be found in the sectors / = 2 and / = 3. We note that while the number of states in 
/ = 1 grows proportionally with L that there are only 4 states with / = 2 and 2 states 



with / = 3. The two / = 2 states of particle 1 at P = 0, 27r/10 are accounted for by 
subtracting ( |2.9| ) from the first equation of (|2.10| ) to write 

oo nC^n-A^'-n nC^n 

(3,4)/ x q 58 -q Bs 

Til,. ..,718=0 

and by considering the ni = 1, n%, . . . ,ng = term of the sum to find (1— g 2 )/ (1— q) = 1+q. 
Similarly, the two / = 2 states of particle 2 at P = 4%/ 10, 67r/10 are obtained from (|3.1|) by 
setting rix = 0, ri2 = 1, 7^3, . . . , n§ = to obtain (q 2 — <? 4 )/(l — q) = q 2 + q 3 - and the two 
I = 3 states of particle 2 at P = 0, 27r/10 are obtained by first using (|2.10|) to write 

00 nC^n-A^-n nC^n-A^'-n 

y (3,4) M _ q E * ~q E * (32) 

77,1 ,...,718=0 

and then setting m = 0, ??,2 = 1, ri3, . . . , rig = to obtain (1 — q 2 )/(l — q) = 1 + q. 

In Fig. 5 we do a similar study for the states which behave as 4A(q) as q — > 1. The 
states which are the single particle excitations for particles 3 and 4 in / = 1 are again 
marked with triangles and of the remaining states there are 9 in / = 2 and 5 in / = 3. As 
before, the character representations (|3J]) and (|3.2| ) account for the / = 3 state of particle 
3 at P = 4tc/10, the / = 2 states of particle 3 at P = 67r/10, 87r/10, 7r, the / = 3 states of 
particle 4 at P = 6tt/10, 8tt/10 and the I = 2 states of particle 4 at P = ir, 12tt/10, 14tt/10. 

We see from Figs. 4 and 5 that the states with / = 2 and / = 3 exactly fill up the 
states missing in I = 1. We also see that the L — > 00 excitation curves give an excellent 
fit to the finite size energies for / = 1 for all momenta which are allowed by ( |2.14|) . Such 
agreement is not particularly seen for the remaining states in I = 2, 3. 

In Figs. 6 and 7 we study this same selection of states for q > and compare with the 
massive dispersion relation ( |2.3| ). Here we see that if q is sufficiently large then all states, 
regardless of their sector, lie on the same massive dispersion curve. 



4. Conclusions 



We may now discuss the findings of |12[-[0 and the relation with [ffl. The study 



in |l^]-[13l is restricted to P = and the major conclusions are 1) that the string content 
of the P = states qualitatively changes with q and 2) that the masses of particles 1 — 5 
as determined from P = agree with the masses of the E$ field theory of @-|§. These 
conclusions are in agreement with our data. On the other hand it is perhaps not exactly 
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fair to say that the single particle states 1 and 3 lie in sector / = 2 and particles 2 and 4 
lie in sector / = 3 because the only sector at q = which contains single particle states 
is I = 1. When q increases there is a motion of the finite number of states in / = 2, 3 to 
smoothly join the macroscopic number of states in I = 1. This is the numerical explanation 
of why the mass computations done in jl| in the sector / = 1 are correct even though P = 
is never in the sector 1 = 1. This motion is seen numerically in ||12||-||13|] where the crossover 



happens on a scale of q ~ L -15 / 8 . It is clearly most desirable to prove this analytically. 

We note that the absence of genuine single particle states in / = 2, 3 means that these 
sectors have a somewhat different nature than I = 1 and it is presumably for this reason 
that a representation of the spectrum in the form (|2.12|) - (|2.16D has not yet been found. 



Such a representation is needed to derive and interpret the character formulas ( |2.10|) . 

We also note that the recognition that the excitations all have an extended Brillouin 



zone scheme renders obsolete the suggestion in [J3[ that the model contains massive 
particles at criticality. 

Finally we remark that the question must be raised as to what the introduction of 
non-integrability (no matter how small) will have on the single particle excitations. In 
particular the pieces of the spectrum in the extended Brillouin zone which lie far up in 
the spectrum and cross a very large number of levels when q —>■ should be expected to 
decay in some fashion when non-integrability is introduced. This should be relevant to 
the relation of the Eg integrable model to the ordinary Ising model in a magnetic field at 
T = T c and remains to be explored. 
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1 3 5 7 8 6 2 
Fig. 1. The Dynkin diagram for the Lie group Eg with the labeling of nodes used in 
this paper. 
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Fig. 2. The eight single particle dispersion relations at qs = 0. 
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Fig. 3. The eight single particle dispersion relations at qs = 0.2. 
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P/tt 

Fig. 4. The plot for qs = and L = 10 of the states which for q — > 1 behave as 
274(5). The states in I = 1 are marked by triangles, in / = 2 by crosses and in / = 3 by 
pentagons. The smooth curves are the theoretical curves ( |2.1| ) for particles 1 and 2 for 
I = 1 and L —>■ 00. Note that particle 1 always lies above particle 2. 
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Fig. 5. The plot for qs = and L = 10 of the states which for q — > 1 behave as 
4A(q). The states in I = 1 are marked by triangles, in / = 2 by crosses and in / = 3 by 
pentagons. The smooth curves are the theoretical curves ( |2.1| ) for particles 3 and 4 for 
I = 1 and L —>■ oo. Note that particle 3 always lies above particle 4. 
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Fig. 6. The plot for qg = 0.2 and L = 10 of the states which for q —> 1 behave as 
2A(q). The states in I = 1 are marked by triangles, in / = 2 by crosses and in / = 3 by 
pentagons. The smooth curves are the theoretical curves ( |2.3| ) for particles 1 and 2 for 
I = 1 and L —>■ oo. Note that these two curves cross near P = n. 
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Fig. 7. The plot for qg = 0.2 and L = 10 of the states which for g — » 1 behave as 
4A(q). The states in I = 1 are marked by triangles, in / = 2 by crosses and in / = 3 by 
pentagons. The smooth curves are the theoretical curves ( |2.3| ) for particles 3 and 4 for 
I = 1 and L —>■ oo. Note that these two curves cross near P = 1.8tt. 
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